Ground state of a resonant two-qubit Rabi model in the ultrastrong coupling regime 
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We consider a generalized Rabi model formed by two identical qubits interacting with a common 
oscillator mode. In the near resonance configuration where the oscillator frequency is close to the 
transition frequency of the qubit, we determine the ground state of the model approximately by 
using a variational method and a transformation method. Both methods are shown to provide good 
agreement with the exact numerical ground state for a range of coupling strength in the ultrastrong 
regime. In addition, we indicate how the accuracy of the approximation can be further improved by 
using perturbation theory. We also examine the ground state entanglement between the two qubits. 
By calculating the negativity, we quantify the degree of entanglement as a function of coupling 
strength. 
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I. INTRODUCTION 

Recently, there has been a revival of interest in study- 
ing the behavior of a two-level system interacting with 
a quantum harmonic oscillator beyond rotating wave ap- 
proximation (RWA). One of the motivations lies in the 
ultrastrong coupling regime recently explored by experi- 
ments with artificial atoms and cavity photon resonators 
1 and microcavities embedding doped quantum wells 
3|, |4j. In such a regime, the qubit-oscillator interacting 
strength is comparable to the oscillator frequency or the 
natural transition frequency of the qubit, and theoretical 
investigations have found novel phenomena, such as the 
asymmetry of vacuum Rabi-splitting km, photon block- 
ade [f| , nonclassical states generation [3] , superradiance 
transition and collapse and revival dynamics 

For a single qubit interacting with a quantum harmonic 
oscillator, the system is described by the Rabi model 
and it has been studied extensively. While there exist an- 
alytical methods [ill Ell to determine the ground state 
exactly, approximation schemes are often employed be- 
cause the approximate ground states in their closed forms 
could be more physically transparent for detailed inves- 
tigations [HHH]. 

In this paper we investigate the ground state of a 
generalized Rabi model in which a quantum harmonic 
oscillator interacts with two identical qubits symmet- 
rically. This is a special case of the Tavis- Cummings 
model [25[ beyond RWA. Recently Agarwal et al. [26| 
have employed an adiabatic approximation to determine 
the energy spectrum and entanglement dynamics in the 
quasidegenerate situation with uja <C w c , where uja and 
oj c are natural frequencies of the qubit and the oscillator 
respectively. Our work here will focus on a different sit- 
uation with uja ~ u c . Such a near resonance situation 
is more complicated because of the break down of the 
adiabatic approximation. 

To approach the problem, we employ two approxima- 
tion schemes to determine the ground state. In Sec. II, we 
describe a variational solution, in which the trial ground 
wave function is constructed from a generalization of a 



single qubit problem [H,[T3|- In Sec. Ill, a transforma- 
tion method is introduced in order to simplify the Hamil- 
tonian. Both methods turn out to give solutions that 
have good agreement with the exact numerical ground 
state. In particular the accuracy of the transformation 
method can be further improved by combining with per- 
turbation theory. 

An interesting feature in the system is the quantum 
entanglement between the two qubits. In Sec. IV, we 
shall characterize the entanglement by using negativity. 
We shall see that there exists entanglement even the sys- 
tem in the ground state. Such an entanglement is due 
to counter-rotating terms involving virtual photons. By 
using our approximate ground state, we quantify the en- 
tanglement by negativity as a function of the coupling 
strength. 



II. THE HAMILTONIAN AND VARIATIONAL 
GROUND STATE 

The Hamiltonian of the model is given by (ft = 1) 

H = loaJx + w c a^ a + g(a + a')J z , (1) 

where a and a) are the annihilation and creation oper- 
ators of the cavity field mode of frequency uj c . The J's 
are spin-1 angular momentum operators. Physically, the 
spin-1 system can be formed by two identical two- level 
atoms in the triplet space. Then uja is the atomic tran- 
sition frequency, and g denotes the collective atom-field 
coupling strength. We note that the (antisymmetric) sin- 
glet state formed by the two atoms does not couple to the 
field and it does not constitute the ground state of the 
system. 

Our task in this paper is to determine the ground state 
energy E g and the ground state vector \ip g ) of the system, 
which are defined by: 



Htyg) = E g \i> g ). 



(2) 



To this end, we let 
quantum number j 



\m) a be an eigenvector of J z (the 
= 1 is suppressed for convenience), 
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FIG. 1: (Color online) [(a)-(b)] Solution of the parameters a and (5 that minimize (H) as a function of g for various detunings. 
(c) The ground state energy as a function of coupling strength g/uiA- The solid line corresponds to E v from the variational 
method and the dashed line corresponds to the exact numerical values E g . (d) The fidelity of the ground state obtained by the 
variational method. 



i.e., J z \m)A — m\rn)A with rn = 0, ±1. Then the 
ground state of the system can be expressed as \i/i g ) = 
Em=o,±i \fm)F\m) A , with \f m ) F being the field state as- 
sociated with the atomic state |m)x- 

To search for |/ m ), we make use of variational method. 
Specifically, we assume that |/ TO ) are coherent states so 
that the trial ground state is in the form: 



\1>v) 



/3|0>f|0>, 



|1>a) 



(3) 

where \a) F is a coherent state of the field with a real 
amplitude a, f3 is a variation parameter, and TV 2 = 2 + /3 2 
is the normalization constant. The expectation value of 
the Hamiltonian (H) with respect to \ip v ) is given by, 



(H) = 



2 + ^ 2 



a 2 oj c 



2ag + V2/3i 



uj A e 



(4) 



By varying the parameters a and /3, (H) is minimized at 
— and = 0, which give the conditions: 



da 



9 



V2 



(auj c - g) e 



aco A 

c?p Uc - ( 2 - /3 2 ) ^ e -4 
2^/3 



(5) 
(6) 



The values of a and /3 can be obtained by solving the 
above equations numerically (Fig. la and b), and hence 
the approximate ground state energy E v — (H) can be 
determined. We remark that there are two sets of a and j3 
which satisfy the conditions given in Eqs. (5) and (6), but 
only the set with a < g/uj c minimizes (H). In addition, 



if the coupling strength g is smaller than co c + uja, then 
a and j3 take the following approximate forms: 



9 



-V2 



2uja + u c 



V2uja (uja + U>c) 



(7) 
(8) 



For the parameter used in Fig. 1, we see that Eqs. (7) 
and (8) can well capture the exact values of a and (3 for 
small g. 

We have tested the accuracy of the variational method 
by comparing the approximate ground state energy E v 
with the numerically exact eigenvalues E g near the reso- 
nance situation with uja ~ w c . The comparison is graph- 
ically shown in Fig. lc. We see that the variational solu- 
tion can produce the ground state energies very well agree 
with the exact values for gjuiA < 0.5. For example, the 
error is about 0.1% at g = 0.5 uja at the exact resonance. 
In fact, the variational solution works in a wider range 
of g when there is a positive detuning uj c — uja > 0. This 
is shown in the curves uj c — 1.2 uja in which there is no 
visible discrepancy for g/ujA < 0.8. 

Another measure of the accuracy is the fidelity defined 
by inner product between the exact numerical ground 
state vector and the variational state vector, i.e., F = 
(tjj v \tpg}. The results are shown in Fig. Id, where we 
have F > 0.999 for g/ujA < 0.5 for the resonance case. 



III. TRANSFORMATION METHOD 

In this section, we describe an alternative method to 
obtain the ground state of the system. Such a method 
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is based on a unitary transformation that simplifies the 
Hamiltonian to a form without countering rotating terms 
approximately. In the case of the (spin-half) Rabi model, 
the corresponding unitary transformation has been dis- 
cussed in Refs. [la - lla ] . Here we indicate how the method 
can be extended to our spin-1 model. 

To begin with we consider a unitary transformation 
defined by H' = e s He~ s . Here 

X (« f ~ a) Jz (9) 



S 



with x being a parameter to be determined. Under the 
unitary transformation, H' can be expressed as 



H> = H> 



H[+H' 2 , 



(10) 



with 



H' = i]uj A J x - (2gx - LOcX 2 ) J 2 + uJc^a, (11) 

H[ = (g - LOcx) (a 1 + a) J z + ir/u A X {a} - a) J y (12) 

H 2 = uj a J x {cosh [x (a f - a)] - r?} 

+ iu A J y {sinh [x (a f - a)] - r\x (a f - a) } • (13) 

Here r/ = e~~ is defined. The cosh \x — a )] an< i 
sinh \x \a) — a)] terms in H 2 has the leading expansion 
term: 

cosh (x (at- a)) =r, + 0{x 2 ), (14) 
sinh (x (a f - a)) = T]X (» f -a)+0 (x 3 ) ■ (15) 

Therefore H 2 is of order 0(x 2 ) describing higher order 
(multi-photon) processes, and we shall neglect H' 2 as an 
approximation, i.e., H' w £Aq + This is the same 
approximation procedure in the spin-half system [16l - ll8j . 
However, the main difference here is the presence of the 
J 2 operator term in H' . This is in contrast to the spin- 
half system in which the corresponding S" 2 term is just a 
constant. 

To proceed, it is useful to express the atomic part of 
Hq+H^ in the diagonal basis of tju)a Jx~ (^9X ~ u cX 2 ) J 2 - 
Such an operator appears in H' a and it describes a 
renormalized three-level system. Specifically, we solve 
[tjujaJx - (2#x - uJcX 2 ) Jl\\v) = e v \v) for eigenvectors 
\v) and eigenvalues e„, with the index v = 0, ± labelling 
the eigenvectors such that e+ > Eq > £-■ It can be 
shown that Eq = —firjWA, e± = tjuja (— A* i \/4 + /x 2 )/2 
with = (2gx — Lo c x 2 ) It}uja, and the eigenvectors are: 



|±} = ^(|i)a + A ± | 



\0) = —(\l) A -\-l) A ), 



a + \-1)a), (16) 
(17) 



where A± = (fj, ± y/4 + n 2 )/y/2 and AT^ = 2 + A| arc the 
normalization constants. In this way, 



H' w w c a) qIa + - 



cia 1 



cia + C2a 
04a + €30 





03a + 04a 



(18) 



where I a is the 3x3 identity matrix operating on the 
atomic subspace. The Cj are some constants depending 
on the parameter x- We see in Eq. (18) that the terms 
involving coefficients c\ and C3 are counter rotating, be- 
cause they describe virtual processes in which atomic ex- 
citation is accompanied by the emission of a photon. 

The main purpose of transforming H into H' is to elim- 
inate some of the counter rotating terms by using a suit- 
able parameter x- Unlike the spin-half atom, it is not 
possible to find a x to remove all counter rotating terms. 
However, since we are interested in the ground state, it 
is sufficient to find x that makes c\ = 0, because in this 
case the corresponding H' in Eq. (18) has the eigenvector 
|0}f|— )a, which is expected to be the ground state vec- 
tor if the interaction is not too strong. Returning to the 
original frame, we have the approximation for the ground 
state energy: 



TjiJA (a 1 + v^ + A* 2 ) 



(19) 



and the ground state vector \tp g ) w \iPt) (subscript T for 
the transformation method) , where 

IVT) = c~ S \0)f\-)a 

1 (lx) F |-l)^ + A-|0) F |Q)^ + |-x)^|l)J. 



(20) 



Here \x)p denotes the coherent state of the field with the 
amplitude x- We see that this is exactly the same form of 
the trial variational function given in Eq. (3). The value 
of x is determined by the condition c\ — which gives, 



V2_ 

X^A 



(X^c - 9) e" 



(21) 



Moreover, by A_ 



^^and. 



2gx " cX , we have 

■qUA ' 



X 2 \-lu c - (2 - A 2 _) 



2xA- 



(22) 



Comparing Eqs. (21) and (22) with Eqs. (5) and (6) ob- 
tained in the variational method, we can see that the 
transformation method and the variational method give 
the same ground state, with x an d A_ are equivalent to 
a and (3 respectively, i.e., we have e_ = E v , \ipr) = \^v)- 
Therefore the transformation procedure in this section 
can be understood as a variation approach, in the sense 
that we vary the Hamiltonian by the unitary transforma- 
tion, instead of the trial wave function. 

However, the main advantage of the transformation 
method is that the corrections can be identified explicitly 
This is because H 2 as is the only approximation made 
in the above derivation. Therefore the exact ground state 
involves higher order processes described by H 2 , and this 
information cannot be obtained by the simply varying 
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the wave function as in Sec. II. Specifically, for x < 1, 
the leading contribution of H' 2 is 



H' 2 



VX 2 ua 



J x (a t2 - 2a 



(23) 
»)] 



which comes from the expansion of the cosh[x (a 
term. By treating Eq. (23) as a weak perturbation, we 
can use the perturbation theory to find the correction of 
the ground state energy (up to second order of H' 2 ): 



E„ 



SE 



e_ + SE, 
N 2 



2ei 



N 2 (2oj c 



(24) 
.(25) 



In Table I, we indicate the performance of the approxi- 
mation scheme for the resonance case uja — ui c . We see 
that although e_ is already quite close to the exact eigen- 
value Eg , there is a significant improvement by including 
the correction 5E, by comparing the second and forth 
columns in the table. 



g/ujA 


Eg/ UJA 


s-/ua 


(e_ + 6E)/ua 


0.2 


-1.01015 


-1.01013 


-1.01015 


0.4 


-1.04256 


-1.04210 


-1.04255 


0.6 


-1.10404 


-1.10137 


-1.10403 


0.8 


-1.20984 


-1.19965 


-1.20988 


1.0 


-1.38986 


-1.36052 


-1.39094 


1.2 


-1.68602 


-1.62699 


-1.68995 



TABLE I: A list of ground state energy obtained by various 



methods at ui c 
in the text. 



UJA- 



Eg, £-, and (e_ + 8E)/uja are defined 



IV. QUANTUM ENTANGLEMENT BETWEEN 
QUBITS 

In this section we examine the quantum entanglement 
between the two qubits when the system is in the ground 
state. Negativity is employed in order to quantify how 
the two qubits are entangled in the ultrastrong coupling 
regime. The negativity is defined as [I?} 



P T A A 



1 



(26) 



where 



Pa 



is the trace norm of p T A , p T A is the partial 

transpose of the reduced density matrix p A = Trj? (p) 
and p = \ifj g ) (ip g \- Alternatively, Eq. (26) can be calcu- 
lated by the absolute value of the sum of the negative 
eigenvalues of p T A A . 

Using the approximate ground state vector \ip v ) in 
Eq. (3), the corresponding reduced density matrix p A 
can be expressed as 



PA ~2N~ 2 



'pn P14 

1 I p 2 2 P23 

P32 P33 

VP41 P44 



:) 



(27) 



where 



Pn = l + /? 2 + 2 V / 2>TTr . 
pu = Pa = 1 - P 2 + e~ 2 ° 2 

P22 = P23 = P32 = P33 = 1 " 

p 44 = l + (3 2 -2V2/3e-^ - 



-2a' 



(28) 
(29) 
(30) 
(31) 



and the order of the columns and rows is |ee), \eg), \ge) 
and \gg) which correspond to the bare atomic level of the 
two qubits. The four eigenvalues of p T A A can be obtained 
analytically and only one of them is negative. Thus, the 
expression of the approximate negativity M is 



M = 



2e 



P 2 



2(2 + /3 2 ) 



(32) 



In the limit g/u>A « 1, a and /3 are approximated by 
Eqs. (7) and (8), and this gives 



TV : 



4uj a (uj a + uj c ) 



(33) 



Therefore the degree of entanglement increases with g 2 
when g/cuA « 1 As a remark, we have also calculated 
the corresponding concurrence and it is simply twice of 
M calculated in the Eq. (32). 

We have computed the negativity from the numeri- 
cally exact ground state for the resonance case, and it 
is plotted in Fig. 2 with a dashed line. We see that the 
negativity is well approximated by Eq. (32) (solid line) 
for g/u)A < 0.5. In particular, the quadratic dependence 
of g is captured by Eq. (33) at small g. There is a greater 
discrepancy when g/uJA > 0.5, which is consistent with 
fidelity F behavior shown in Fig. Id. Although Fig. 2 
does not show the case of g/uJA > 1, the exact negativity 
drops as g /uja further increases after it reaches its maxi- 
mum at g ss uja- When g rj 2.6 uja, the negativity would 
drop to zero and never increase again. In other words, 
the two qubits become disentangled if g is sufficiently 
large. 
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J\f 0.02 



0.01 
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FIG. 2: (Color online) The negativity j\f as a function of 
coupling strength g/u) A for the ui c = uja case. The solid line 
(dashed line) corresponds to the negativity of the approximate 
(exact) ground state. 
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V. CONCLUSION 

To conclude, we have shown that the ground state of 
the Rabi model formed by two identical qubits in the 
near resonance configuration can be determined approx- 
imately by using a variational method and a transfor- 
mation method, and the results are in good agreement 
with exact numerical calculations when g is a significant 
fraction of loa (up to g/u)A = 0.5). For instance, the er- 
ror of the variational ground state energy for the exact 
resonance case is about 0.1% at g = 0.5 uja, and the er- 
ror can be greatly reduced to 0.0004% by further using 
second order perturbation theory. The key advantage of 
our methods is that the analytical form of the ground 
state can be captured approximately, which is useful for 
analyzing the ground state properties. Specifically, we 



have examined the quantum entanglement between the 
two qubits. Such an entanglement is a consequence of 
counter rotating terms in the Hamiltonian, because if 
such terms are dropped by RWA, then the two qubits 
are simply disentangled in the ground state. Using our 
approximate ground state vector, we are able to deter- 
mine the negativity as a function of system parameters 
for g/uiA up to 0.5. 
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